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Abstract 

We calculate the flux threading the pick-up coil of a square SQUID magnetometer 
in the presence of a current dipole source. The result reproduces that of a circle coil 
magnetometer calculated by Wikswo [1] with only small differences. However it has 
a simpler form so that it is possible to derive from it closed form expressions for the 
current dipole sensitivity and the spatial resolution. The results are useful to assess 
the overall performance of the device and to compare different designs. 
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1 Introduction 



The application of Superconducting QUantum Interference Devices (SQUIDs) 
to the measurement of biomagnetic fields has occurred because of their sensi- 
tivity, their stability and their flexibility. In fact SQUID based devices offer the 
Email address: s.rombetto@cib.na.cnr.it (S Rombetto ). 



Preprint submitted to Elsevier 



31 March 2008 



possibility to implement measurements where no other methodology is possi- 
ble and moreover they present the advantage to be a non-invasive technique 
[2J, [3]. Usually studies, aimed to develop better devices, focus on gradiometric 
configurations since they are less sensitive to the noise. Hereafter we focus on 
magnetometer configuration which is successfully employed in multichannel 
systems for biomagnetic imaging [3]. Magnetometers are extremely sensitive 
to the outside environment, while some other configuration, like gradiometers 
provide the advantage of discriminating against unwanted background fields 
from distant sources while retaining sensitivity to the nearby sources. 
In a dc-SQUID magnetometer, the pick-up coil, collects the magnetic flux giv- 
ing an effective area much larger than that of the SQUID itself [5] . Here we 
study the effect of magnetometer pick-up coil geometry on the performances 
of SQUID devices for biomagnetism. Well known and widely spread expres- 
sions for the flux threading a magnetometer and the current dipole sensitivity 
have been calculated by Wikswo [1] referring to a circular loop device. Since 
typically, SQUID magnetometers present a square pick-up loop [5], we were 
driven, for the best characterization of such devices, but also for general rea- 
sons, to recalculate the quantities of interest in the case of square pick-up loop. 
It turned out that expressions for the minimum detectable current dipole and 
for the spatial resolution are easily derived for the square geometry. 

2 Magnetic flux threading a square magnetometer in the current 
dipole model 

A widely used mathematical model to describe bioelectric currents is the cur- 
rent dipole [6] . It is a good model of elementary cellular events, thus it can be 
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used for magnetoencephalography as well for magnetocardiography studies. 
Let us consider a dipolar electric current source p(p x ,P y ) located at a point 
r' = (x', 0, z') in a conducting half space and a pick-up loop centered on 
the z-axis (Fig{TJ). The magnetic field generated by the source p at a point 
r = (x, y, z) has the vector potential A (A x , A y , 0) 




Fig. 1. Schematic of the square loop magnetometer. A current dipole source p is 
placed at the point r' = (x' , 0, z'). 

By using eq.([T]) it is possible to calculate the magnetic flux through the con- 
sidered pick-up coil by performing a line integral around the loop. In the case 
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of a circular loop having radius R, Wikswo pQ derived the result 
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K(k) and E(k) are respectively the complete elliptic integrals of first and 
second kind. In a very similar way we have calculated the flux collected by a 
square loop magnetometer having size L, laying parallel to the xy plane, at a 
distance D = z — z' above a current dipole source p (FigJT]). The result is 
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Note that p x does not contribute to the collected flux for symmetry rea- 
sons, so that there is no loss of generality if we consider p as having only 
the y— component. In order to compare devices presenting different geome- 
tries (square loop or circle loop), we shall consider equal area devices, that 
is equivalent to the condition R = L/^/tt. Also we shall introduce in the 
above equations the dimensionless source position x = x' / D, the reduced flux 
^7i/fi p y and the geometrical parameter q = L/D. As we shall see, although 
eq.pj) and eq.(jl]) give, almost, the same result, eq.(pl| allows for useful analyt- 
ical progresses which eq.(E]) does not permit. 
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3 Maximum magnetic flux 



In order to evaluate the minimum detectable current dipole and the spatial 
resolution as a function of the geometrical parameters L and D, it is essential 
to determine the value of the maximum flux and the maximizing source posi- 
tion x' = x' max for any q. 

In FigfS]the reduced flux, calculated by means of eqs.fl2]), (jl]), is plotted as a 
function of the reduced source position x for two different values of the ratio 
q, for the two geometries, circle and square. As a general picture one sees that 
the flux is zero when the source is exactly under the loop (V = 0) and, as the 
source moves away, it maximizes for x' = x' max , before decaying out. 
While the elementary procedure (zeros of the derivatives) for the maximum 
finding does not give straightforward results for eq.([5]) and eq.fllj), simple ap- 
proximate analytical results can be obtained for the maximum flux, on the 
basis of physical considerations. First consider the case of a large loop size to 
source distance ratio (q 3> 1). It is evident that in this case (in exact manner 
in the limit of zero distance) the collected flux is maximum when the source 
position coincides with the loop edge. This is to say x' = L/2 (x = q/2) and 
D — > for a square loop, and analogously, x' = R (x = q/\pn) and D — > 
for a circle loop. Thus we can take this asymptotic values, q/2 and q/y/ir, 
as approximate maximum positions for the circle and the square loop, as it 
is shown in Figj2j, where the two maxima correspond roughly to the points 
x = 2.5 and x = 2.82 respectively, since q=5 for the curves in Fig.(T5]). 
In the case of a square loop, using the value q/2 in eq.(j3J, we are lead to the 
following result for the maximum flux 
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Now we turn to the case when the ratio between the loop size and the distance 
is quite small (g < 1). Equations ([2]) and (j4j) give in practice identical results 
for the two loop geometries and the difference between the fluxes threading 
devices with different shape cannot be appreciated: the two results overlap. 
For large distances above the source, or small loop area, equations ([S]),® can 
be simplified by substituting R = L/^fn in eq.fl2J), and then expanding these 
expressions in the small parameter q. Both equations give the same result 

$ = VoPy x /L\ 2 _ jJ, Py X (R 
4 (l + x 2 f 2 \D) ~ 4vr (l + x 2 f 2 \D 

which is shown in Figj2]by the gray curves on the left, calculated for q = 0.5. 
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Fig. 2. Comparison between the normalized flux threading devices presenting square 
or circle shape as a function of the source position x'/D. Gray curves on the left 
have been obtained using q = 0.5 in eq. The two black curves on the right have 
been obtained using q = 5 in eq. The gray dotted curve shown on the left has 
been obtained by using for both geometries the Taylor expansion, given by eq. 
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thus eq.Q represents the flux as a function of x, threading a magnetometer 
positioned at large distance, with respect to the loop size, above the source. 
In FigfSJthe dotted curve represents results obtained by eq.flHD- In this regime 
(q <C 1) circle and square loop, give identical results, so that any information 
about the shape of the magnetometer loop is lost. It is easily found that eq.© 
maximizes exactly for x — (x' — ^D) so that the maximum flux, for small 
loop size to source distance ratio (q <C 1) is 



Thus small differences in the collected flux, due to the inhomogeneity of the 
source, emerge between square and circle geometry only when the source is 
close faced to the loop. This situation is illustrated in Figf5]by the black curves 
on the right, calculated for q = 5. As can be seen, the two curves maximize in 
slightly different points, as already observed. 

A quantitative evaluation, based on numerics, of the validity of the approxi- 
mations given by eq.([5]) and eq.([7|) as well as of the maximizing positions will 
be given in the next section. We close this section by an estimation of the 
maximum flux for a typical situation. If we consider for the current dipole the 
value p y = WnA ■ m we find that $^ re = 9 -18 • 10 -17 Wb = 45 m$ for 
q = 0.5 (obtained by using the values L = 9 mm, D — 1.8 cm). 

4 Minimum detectable current dipole 

In order to determine the smallest detectable current dipole p r y nm , we have to 
impose that the collected flux $^™ re is comparable to the total flux noise <3>*. 
Therefore p™ n is an evaluation for the sensitivity of the considered device: if 
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a device can detect a smaller p y , then it presents a higher sensitivity. 
For large q values (q ^> 1), from eq. (JSJ) for x = q/2 we obtain 
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For small g, from eq.© we obtain 

f = ^ (9) 



The dependence of the sensitivity on the loop to source distance ratio de- 
scribed by eqs.flSD, @ is shown in FigJHJ In the same figure it is also shown 
for a comparison, p™ m evaluated by a numerical maximum finding procedure 
directly from eq.flU). 

For very small loop size to source distance ratio (g— > 0) the current dipole 
sensitivity diverges as q~ 2 , due to the small area of the SQUID pick-up loop. 
In the opposite limit, i.e. for very small distance between source and sensor or 
very large SQUID sensors (q— > oo), the sensitivity improves without limits 
{py— > oo) because the collected flux continues to grow. 



5 Spatial resolution 

When the current dipole source moves from the position of the maximum flux 
along the x direction, with a displacement 5, there is a change in the flux A$. 
In a general way, for small 5, one obtains the following expression 

A$ = $( x ')| x , +s - <f>(x')L, = 
H Li 5 2 )-$(x%, = — L 5 2 (10) 



8 



where x' max is the value for which the flux maximizes and the condition 
&(x' max ) = has been used. 

If we now assume that the spatial resolution is "the least detectable displace- 
ment" corresponding to a variation in flux equal to the flux noise $*, by 
inverting eq. ffTOj) one obtains 
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The smaller is S, the better is the resolution. 

We now derive an analytical expression for the spatial resolution. For large 




q=L/D 



Fig. 3. The spatial resolution and the current dipole sensitivity versus the ratio of 
the side of a square sensor to the source-sensor distance. Black curves describe the 
current dipole sensitivity: solid curve is for the case q S> 1 (eq\Bjj, dotted curve is 
for the case q <C 1 ( egHjj and circles are for the numerical calculation. Gray curves 
describe the spatial resolution: solid curve is for the case g > 1 (eq \12\i , dashed 
curve is for the case q <C 1 ( eq [T3\) and squares are for the numerical calculation. 
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loop size to source distance ratio (q ^> 1), as we have seen before, the flux 
maximizes approximatively when the condition x' max (q) = q/2 is satisfied. In 
this regime, the expression derived for spatial resolution from eqs.(jl]), fllip . 
developed around x — | , becomes 
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For small loop to source distance ratio (q <C 1), the analytical expression for 
the spatial resolution can be obtained on the basis of eqs.fjSJ) and (TTTT) . for 
x = and the expression for spatial resolution is 
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When the g value is about 1, the approximations introduced till for 6 now begin 
to fail, so that it is necessary to compute the spatial resolution numerically. 
In order to do this and for a comparison with the analytical results, we have 
calculated analytically the second derivative of the flux given in eq.(jH) and 
calculated its value in the maximizing position x' max evaluated numerically for 
any q. 

In Figj3] the two analytical solutions eq. (fT2|) and eq. (fT3]) .and the numerical 
result for the spatial resolution are plotted in gray. Solid curve is for the case 
q 1, dashed curve is for the case gCl and squares are for the numerical 
calculation. 

It is worth noting that the spatial resolution 5 defined in eq. ffl2|) has the lower 
limit (obtained for q — > oo) 
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This means that even if we design a device that could collect a very large flux, 
the spatial resolution cannot enhance. 

6 Conclusions 

We found expressions for the dipole sensitivity and for the spatial resolution of 
a square loop magnetometer starting from eq.(j4"j) which gives the flux thread- 
ing the loop, due a current dipole source. Both quantities show a monotonic 
dependence on the sensor size for a fixed sensor to source distance. The dipole 
sensitivity is limited only by the loop size L. On the contrary the spatial res- 
olution has a lower limit given by eq. (|14p . meaning that there is no way to 
improve the spatial resolution even using a very large loop size device. Thus for 
all practical needs the calculations here presented indicate that when the dis- 
tance D is comparable with the size of the loop L, the limit spatial resolution 
is already obtained. 
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